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Conditional  Sampling  in  a 
Transitional  Boundary  Layer 
Under  High  Freestream 
Turbulence  Conditions 

Conditional  sampling  has  been  performed  on  data  from  a  transitional  boundary  layer 
subject  to  high  (initially  9%)  freestream  turbulence  and  strong  (K  =  (v/Ul,)(dU„/dx)  as 
high  as  9X10~6)  acceleration.  Methods  for  separating  the  turbulent  and  nonturbulent 
zone  data  based  on  the  instantaneous  streamwise  velocity  and  the  turbulent  shear  stress 
were  tested  and  found  to  agree.  Mean  velocity  profiles  were  clearly  different  in  the  tur¬ 
bulent  and  nonturbulent  zones,  and  skin  friction  coefficients  were  as  much  as  70%  higher 
in  the  turbulent  zone.  The  streamwise  fluctuating  velocity,  in  contrast,  was  only  about 
10%  higher  in  the  turbulent  zone.  Turbulent  shear  stress  differed  by  an  order  of  magni¬ 
tude,  and  eddy  viscosity  was  three  to  four  times  higher  in  the  turbulent  zone.  Eddy 
transport  in  the  nonturbulent  zone  was  still  significant,  however,  and  the  nonturbulent 
zone  did  not  behave  like  a  laminar  boundary  layer.  Within  each  of  the  two  zones  there  was 
considerable  self-similarity  from  the  beginning  to  the  end  of  transition.  This  may  prove 
useful  for  future  modeling  efforts.  [DOI:  10. 1 1 15/1. 1521957] 


Introduction 

Boundary  layer  transition  is  an  important  phenomenon  experi¬ 
enced  by  the  flow  through  gas  turbine  engines.  Mayle  [1]  stated 
that  a  substantial  fraction  of  the  boundary  layer  on  both  sides  of  a 
gas  turbine  airfoil  may  be  transitional.  The  extended  transition 
zones  exist  due  to  strong  favorable  pressure  gradients,  found  on 
both  the  pressure  side  and  the  leading  section  of  the  suction  side, 
which  stabilize  the  boundary  layer  and  delay  transition  in  spite  of 
the  high  freestream  turbulence  intensity  (FSTI)  in  gas  turbine  en¬ 
vironments.  The  ability  to  model  and  predict  high  FSTI  transition 
is  important  since  heat  transfer  rates  and  skin  friction  coefficients 
may  increase  substantially  when  a  boundary  layer  undergoes  tran¬ 
sition.  Boundary  layer  separation,  which  is  believed  to  be  a  sig¬ 
nificant  problem  on  the  suction  side  of  some  low-pressure  turbine 
airfoils,  also  depends  strongly  on  the  state  of  the  boundary  layer 
with  respect  to  transition.  Improved  transition  models  and  turbine 
designs  depend,  therefore,  on  a  better  understanding  of  high  FSTI 
transition. 

Documentation  of  high  FSTI  transition  has  included  work  in 
turbine  cascades  and  rotating  rigs  (e.g.,  Halstead  et  al.  [2]).  Sev¬ 
eral  studies  have  considered  the  flow  along  flat  and  curved  walls. 
These  simpler  geometries  allow  for  more  detailed  in-flow  mea¬ 
surements  than  are  typical  for  cascade  studies.  On  surfaces  subject 
to  zero  streamwise  pressure  gradients,  Blair  [3],  Sohn  and  Re- 
shotko  [4],  and  Kim  et  al.  [5]  all  showed  that  at  FSTI  above  about 
3%,  transition  occurred  rapidly  near  the  leading  edge  of  a  test 
surface.  Blair  [6]  considered  cases  with  FSTI  up  to  5%  and  con¬ 
current  acceleration,  holding  the  acceleration  parameter,  K,  con¬ 
stant  at  values  up  to  0.75X1CP6.  The  acceleration  delayed  the 
transition,  even  with  high  FSTI.  Volino  and  Simon  [7-9]  consid¬ 
ered  transition  along  a  concave  wall  with  inlet  FSTI  of  8%  and 
acceleration  with  K  as  high  as  9X  10~6.  Acceleration  rates,  Rey¬ 
nolds  numbers  and  FSTI  were  typical  of  the  pressure  side  of  a  gas 
turbine  airfoil.  An  extended  transition  region,  with  intermittent 
turbulent  and  nonturbulent  zones,  covered  most  of  the  test  surface. 
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The  turbulent  zone  included  a  range  of  both  large  and  small  scale 
eddies,  much  like  a  fully  turbulent  boundary  layer  and  the  turbu¬ 
lent  zone  in  a  low  FSTI  transitional  boundary  layer.  The  nontur¬ 
bulent  zone  and  the  pre-transitional  boundary  layer  were  not 
laminar-like  as  in  a  low  FSTI  flow,  but  instead  were  characterized 
by  high-amplitude  large-scale  fluctuations  and  an  absence  of 
smaller  scales.  Nonturbulent  zone  velocity  fluctuations  are  be¬ 
lieved  to  be  induced  by  the  freestream  unsteadiness  through  pres¬ 
sure  fluctuations,  as  discussed  in  Volino  [10].  Near-wall  turbu¬ 
lence  production  is  believed  to  be  largely  absent  in  the 
nonturbulent  zone. 

The  intermittent  nature  of  transition,  both  at  low  and  high  FSTI, 
has  led  to  efforts  to  incorporate  intermittency  in  transition  models 
and  to  model  the  two  zones  of  the  intermittent  flow  separately. 
Among  recent  efforts  are  the  work  of  Steelant  and  Dick  [11], 
Suzen  and  Huang  [12],  and  Solomon  et  al.  [13].  Separate  model¬ 
ing  of  the  nonturbulent  and  turbulent  zones  requires  knowledge  of 
the  flow  behavior  within  each  zone,  which  can  be  provided 
through  conditional  sampling  of  experimental  data.  Conditional 
sampling  results  have  been  presented  for  zero  streamwise  pressure 
gradient  conditions  by  Kim  et  al.  [5],  Sohn  and  Reshotko  [4],  and 
Kuan  and  Wang  [14].  Results  from  favorable  pressure  gradient 
cases  have  been  reported  by  Blair  [6],  Wang  and  Keller  [15],  and 
Wang  and  Zhou  [16].  The  strongest  acceleration  rate  in  these 
cases  was  K  =  0.75X  10~6.  This  is  a  relatively  mild  acceleration 
for  modern  airfoils,  and  the  transition  zones  were  short  in  all  cases 
with  elevated  FSTI. 

Volino  and  Simon  [7,9]  obtained  rough  estimates  of  the  time- 
averaged  intermittency  (fraction  of  time  the  flow  is  turbulent)  in 
their  high  FSTI,  strong  acceleration  case  using  an  analog  circuit, 
which  determined  when  the  flow  was  turbulent  based  on  the  time 
derivatives  of  hot-wire  voltages.  The  circuit  worked  well  in  the 
low  FSTI  cases  of  Kim  et  al.  [5],  but  in  high  FSTI  cases  differ¬ 
ences  in  fluctuation  level  between  the  turbulent  and  nonturbulent 
zones  are  narrower,  and  more  careful  setting  of  thresholds  and 
smoothing  of  the  intermittency  function  are  required  for  condi¬ 
tional  sampling.  This  is  difficult  with  an  analog  circuit,  and  is 
better  done  in  post-processing  after  the  signal  has  been  digitized. 
Such  post-processing  requires  data  acquisition  at  a  high  enough 
sampling  rate  to  provide  an  essentially  continuous  signal.  Due  to 
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data  storage  limitations  at  the  time,  Volino  and  Simon  [7,9]  only 
sampled  at  a  high  enough  rate  for  such  processing  at  a  few  iso¬ 
lated  points  in  the  boundary  layer.  Volino  [17]  examined  the  data 
from  these  points  using  wavelet  analysis  and  presented  prelimi¬ 
nary  results  of  conditional  sampling.  Uncertainties  were  high  due 
to  the  limited  amount  of  data.  To  the  authors’  knowledge,  no  other 
detailed  conditional  sampling  results  from  boundary  layers  with 
strong  acceleration,  extended  transition  zones  and  FSTI  greater 
than  5%  are  available  in  the  literature. 

The  present  paper  has  two  objectives.  The  first  is  to  provide  the 
results  of  conditional  sampling  for  a  transitional  boundary  layer 
with  high  FSTI  and  strong  acceleration.  The  case  presented  by 
Volino  and  Simon  [7-9]  has  been  reproduced  on  a  flat  test  wall, 
and  data  have  been  acquired  throughout  the  boundary  layer  at  a 
sufficiently  high  sampling  rate  for  conditional  sampling  post  pro¬ 
cessing.  The  second  objective  is  to  provide  a  baseline  case  for  a 
study  of  the  significance  of  streamwise  curvature  on  transitional 
boundary  layers  at  elevated  FSTI.  A  second  paper  (Schultz  and 
Volino  [18])  presents  results  from  an  otherwise  similar  case  on  a 
wall  with  strong  concave  curvature. 

Experiments 


primarily  due  to  the  increasing  freestream  velocity.  The  stream- 
wise  component  u'x  does  decay  somewhat  due  to  straining  of  the 
freestream  eddies  in  the  accelerating  flow,  but  v  remains  nearly 
constant  at  all  stations. 

Velocity  profile  measurements  were  made  at  ten  streamwise 
stations  along  the  centerline  of  the  test  wall  using  a  single  sensor, 


Table  1  Flow  parameters  at  measurement  stations 
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Facility  and  Measurements.  All  experiments  were  con¬ 
ducted  in  the  low  speed  wind  tunnel  shown  in  Fig.  1.  Three  fans 
supply  air  to  a  plenum,  which  is  followed  by  a  diffuser,  a  settling 
chamber  containing  a  honeycomb,  a  screen  pack,  a  second  settling 
chamber,  and  a  three-dimensional  contraction  which  reduces  the 
cross  sectional  area  from  1.09  mX1.09  m  to  0.69  mX0.18  m.  A 
bi-planar  turbulence-generating  grid  with  43%  blockage  is  in  the 
plane  of  the  contraction  exit.  The  grid  is  based  on  the  design  used 
by  Kim  et  al.  [5]  and  consists  of  a  3.8-cm  diameter,  0.69-m  long 
vertical  pipe  down  the  center  of  the  contraction  exit,  and  six 
3.2-cm  diameter,  0.1 8-m  long  evenly  spaced  horizontal  pipes.  Fol¬ 
lowing  the  grid  is  a  1-m  long  rectangular  development  section  and 
the  test  section,  a  converging  channel.  One  side  of  this  channel  is 
a  flat  Plexiglas  plate  of  0.69  m  width  and  1.2  m  length,  which 
serves  as  the  test  wall.  Pressure  taps  are  installed  along  its  span- 
wise  centerline.  At  the  leading  edge  of  the  test  wall  a  slot  is  used 
to  bleed  off  the  boundary  layer  which  grows  in  the  development 
section.  Opposite  the  test  wall  is  a  flexible  wall  which  can  be 
adjusted  to  set  the  desired  pressure  gradient  along  the  test  wall. 
For  the  present  study  the  inlet  velocity  is  set  to  4.6  m/s  and  the 
velocity  gradient  along  the  wall  is  held  constant  at  13.9  s_I.  The 
acceleration  parameter,  K ,  drops  from  a  maximum  of  9X  10“  6  at 
the  inlet  to  the  test  section  to  1XKV6  at  the  last  measurement 
station.  Values  of  K  through  the  test  section,  measurement  loca¬ 
tions  and  other  parameters  are  given  in  Table  1 . 

At  the  inlet  to  the  test  section  the  mean  flow  is  spatially  uniform 
to  within  3%  and  the  turbulence  is  uniform  to  within  6%.  The 
components  of  the  freestream  turbulence  intensity  are  8.8%,  8.9%, 
and  8.3%  in  the  streamwise,  cross-stream,  and  spanwise  direc¬ 
tions,  respectively.  The  integral  length  scales  of  these  components 
of  the  freestream  turbulence  are  3  cm,  1.6  cm,  and  1.4  cm.  Spectra 
of  the  freestream  turbulence  at  the  most  upstream  measurement 
station  are  shown  in  Fig.  2.  In  the  test  section,  the  freestream 
turbulence  intensity  (normalized  using  the  local  freestream  veloc¬ 
ity)  drops  to  about  2%  at  the  last  measurement  station.  The  drop  is 
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Fig.  2  Freestream  spectra  at  Station  1 
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boundary  layer  type  hot-wire  probe  (TSI  model  1218-T1.5),  a 
boundary  layer  cross-wire  probe  (TSI  model  1243-T1.5),  and  a 
constant-temperature  hot-wire  anemometer  (TSI  model  IFA-100). 
The  probes  were  moved  normal  to  the  wall  at  each  station  using  a 
motorized  traverse  with  minimum  step  size  of  12.5  gm.  At  each 
position  in  the  velocity  profiles,  data  were  acquired  for  26  seconds 
at  a  20  kHz  sampling  rate  (2 19  data  points).  The  hot-wire  signals 
were  low-pass  filtered  at  10  kHz.  All  raw  data  were  stored.  Sam¬ 
pling  at  20  kHz  provided  essentially  continuous  velocity  traces  for 
subsequent  processing. 

Data  Processing.  Mean  and  rms  fluctuating  velocities  in  the 
streamwise  and  wall-normal  directions  were  computed  from  the 
instantaneous  data.  Uncertainties  in  these  quantities  are  3-5%  ex¬ 
cept  in  the  very  near  wall  region  (y  +  <5)  where  near-wall  cor¬ 
rections  (Wills  [19])  are  applied  to  the  mean  velocity.  Uncertainty 
in  the  turbulent  shear  stress,  —u'v',  is  10%.  Skin  friction  coeffi¬ 
cients  were  determined  using  a  technique  whereby  tk  was  ad¬ 
justed  until  the  U+  versus  y+  data  fit  profiles  computed  using 
near-wall  similarity,  as  described  by  Volino  and  Simon  [20].  This 
technique  accounts  for  pressure  gradient  effects  on  the  near  wall 
profile.  Uncertainty  in  Cf  is  8%.  Boundary  layer  thicknesses  were 
determined  from  the  mean  velocity  profiles.  Uncertainties  in  the 
momentum  and  displacement  thicknesses  are  10%.  These  uncer¬ 
tainties  include  bias  uncertainties  which  tend  to  cancel  such  that 
the  uncertainty  in  the  shape  factor,  H,  is  7%. 

Intermittency  Based  on  u.  The  intermittency  function,  T,  in¬ 
dicates  whether  the  boundary  layer  is  instantaneously  turbulent  or 
nonturbulent  at  a  measurement  location.  It  is  assigned  a  value  of 
zero  for  nonturbulent  flow  and  one  for  turbulent  flow.  The  time 
average  of  T  is  the  intermittency,  y.  Keller  and  Wang  [21]  and 
Solomon  [22]  review  several  techniques  for  determining  T.  In  the 
present  study,  two  intermittency  detection  techniques  were  uti¬ 
lized  and  compared.  The  first,  as  used  in  Volino  and  Hultgren  [23], 
is  based  on  the  instantaneous  streamwise  velocity,  u.  The  u  signal 
is  first  digitally  high-pass  filtered  with  a  cutoff  frequency 

fHP  =  Wto-U„  (1) 

where  fHP  is  in  Hz  and  U „  is  in  m/s.  The  filter  eliminates  low- 
frequency  fluctuations,  which  are  common  to  both  the  turbulent 
and  nonturbulent  zones.  The  filter  was  used  by  Volino  and  Hult¬ 
gren  [23]  to  remove  narrow  band  fluctuations  associated  with 
shear  layer  instabilities  in  separated  boundary  layers.  In  the 
present  work,  this  filter  had  little  effect  on  the  intermittency  func¬ 
tion.  The  filtered  signal  is  then  used  to  determine  an  intermittency 
function  Tlu ,  as 

1  if  \du/dt\>i.%9-Um-Um 
r  it)  —  m  n't 

lu3  '  0  otherwise  (  2 

where 

U  +  2-^1/3-7Il/2  if  2-^7Tl/3>^7Il/2 
Um=  .  (3) 

U  otherwise 

The  velocity  U  is  the  local  mean  velocity  at  the  measurement 
location  and  is  the  local  freestream  velocity  at  the  measure¬ 
ment  station.  All  velocities  are  in  m/s,  and  the  time  derivative  is  in 
m/s2.  The  time  derivative  of  the  velocity  should  scale  with  the 
magnitude  of  the  velocity  and  the  frequency  with  which  eddies 
pass,  which  in  turn  also  scales  with  the  velocity.  This  explains  the 
velocity  squared  term  in  the  threshold  in  Eq.  (2).  The  use  of  Um 
instead  of  £/„  allowed  the  threshold  level  to  adjust  in  the  near¬ 
wall  region  where  the  mean  velocity  becomes  small.  Very  near  the 
wall,  as  U  approaches  zero,  the  threshold  would  also  approach 
zero  and  T  l  u  would  go  to  1  if  Um  were  set  equal  to  U.  To  prevent 
this,  Um  is  adjusted  using  the  fluctuating  velocity,  as  shown  in  Eq. 
(3).  This  adjustment  models  the  instantaneous  rise  in  the  near  wall 
velocity  which  occurs  during  turbulent  intervals  when  higher 
speed  fluid  sweeps  toward  the  wall. 


A  second  function,  T2„ ,  is  next  computed  based  on  the  abso¬ 
lute  value  of  the  second  derivative  of  the  filtered  velocity  signal, 
\d2u/St2\  .  The  threshold  for  T  2„  is  set  such  that  the  time  average 
of  r,„  and  T2„  are  equal.  Next,  a  combined  intermittency  func¬ 
tion,  T 3u  is  defined  as 


1  ifri=lorr2=l 
0  otherwise 


(4) 


Basing  T3„  on  both  the  first  and  second  time  derivative  of  the 
velocity  helps  to  minimize  dropouts  within  the  turbulent  zone, 
which  occur  in  Tlu  and  T 2„  when  the  derivatives  cross  zero.  The 
final  step  is  to  smooth  T3„  to  minimize  false  turbulent  points  in 
the  nonturbulent  zone  and  false  nonturbulent  points  in  the  turbu¬ 
lent  zone.  The  function  T3„  is  low-pass  filtered  with  cutoff 
frequency 


fLP=  17.78-  U,  (5) 

where  fLP  is  in  Hz  and  t/„  is  in  m/s.  Finally,  using  the  filtered 
r3„,  the  intermittency  function  T„  is  set  as 


r,(0 


i 

o 


if  r,„>o.5 

otherwise 


(6) 


The  thresholds  and  filter  frequencies  presented  above  were  set 
after  visual  inspection  of  many  data  traces  and  the  T„(r)  resulting 
with  several  different  thresholds.  Visual  inspection  of  the  data  is 
ultimately  the  best  criteria  available  for  determining  how  well  an 
intermittency  function  is  separating  the  turbulent  and  nonturbulent 
flow  [17,22],  The  thresholds  and  frequencies  were  useful  both  for 
the  present  study  and  the  separated  flow  transition  study  of  Volino 
and  Hultgren  [23].  They  are  not  expected  to  be  universal  criteria 
for  all  flows,  and  other  intermittency  detection  schemes  might 
work  as  well.  The  thresholds  might  be  made  more  general  if  non- 
dimensionalized  using  the  viscosity  or  a  characteristic  length  such 
as  the  boundary  layer  thickness. 


Intermittency  Based  on  u'v’  ■  A  second  intermittency  func¬ 
tion  is  computed  based  on  the  instantaneous  turbulent  shear  stress. 
Without  pre-filtering  of  the  velocity  signal,  an  intermittency  func¬ 
tion  is  computed  as 


1 

ri«B(0  =  0 


if  (du'v'ldt)2>50-  Ui 
otherwise 


(V) 


A  second  function  T 2uv  is  computed  based  on  {d2u'  v'  /  St2)2  with 
the  threshold  set  such  that  the  time  averages  of  Tlui,  and  T 2uv  are 


Fig.  3  Instantaneous  velocity  traces  and  intermittency  func¬ 
tion:  (a)  uand  r„,  (b)  u'v'  and  l'„„ 
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Fig.  4  Intermittency  profiles  based  on  (a)  u  and  (to)  Reynolds 
shear  stress 


Fig.  5  Peak  intermittency  in  profile  versus  dimensionless 
streamwise  location 


equal.  The  remaining  steps  in  computing  the  intermittency  func¬ 
tion  r„„(f)  follow  the  steps  in  Eqs.  (4)-(6)  above. 

Figure  3  presents  a  typical  signal  from  the  transitional  boundary 
layer  along  with  both  intermittency  functions.  Some  regions  are 
clearly  turbulent  (e.g.,  0.025  s-0.035  s)  and  some  are  clearly  non- 
turbulent  (0.01  s-0.02  s).  Between  zones  the  demarcation  is  not 
always  sharp  (0.04  s-0.045  s)  leading  to  differences  in  T„(t)  and 
T„„(f)  and  the  possibility  of  some  “leakage”  of  data  between 
zones.  The  correlation  coefficient  between  T„(f)  and  T„„(t)  was 
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always  at  least  0.8,  however,  and  at  95%  of  the  measurement 
locations  it  was  above  0.9.  Although  the  shear  stress  is  often  con¬ 
sidered  a  better  criterion  function,  the  good  agreement  between 
T„(0  and  T„„(0  suggest  that  both  are  acceptable.  This  is  particu¬ 
larly  useful  very  near  the  wall,  where  only  u  can  be  measured. 
Both  functions  were  used  in  conditional  sampling,  providing  es¬ 
sentially  the  same  results.  The  uncertainty  in  the  intermittency,  y, 
is  0.1.  In  the  figures  which  follow,  data  points  are  shown  for  the 
nonturbulent  zone  only  when  the  local  y<95%,  and  for  the  tur- 


25 


1  10  100  1000 


y* 


Fig.  6  Mean  velocity  profiles  in  wall  coordinates:  (a)  compos¬ 
ite;  (to)  nonturbulent;  (c)  turbulent 
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bulent  zone  only  when  y>5%.  Bulk  parameters  such  as  Cy  are 
presented  when  ypk<95%  and  ypk>5%  for  the  nonturbulent  and 
turbulent  zones,  respectively. 


Results 


Intermittency.  Intermittency  profiles  for  the  ten  measure¬ 
ment  stations  are  shown  in  Fig.  4.  Agreement  between  the  u  and 
shear  stress  based  y  is  good.  Intermittency  remains  low  for  the 
first  three  stations  and  then  begins  to  rise.  By  the  last  station, 
transition  is  nearing  completion.  Transition  criteria  such  as  those 
presented  by  Johnson  [24]  and  Mayle  [1]  indicate  that  the  present 
boundary  layer  would  be  transitional  by  Station  1  under  zero  pres¬ 
sure  gradient  conditions.  In  fact,  the  intermittency  is  nonzero  at 
Station  1,  but  the  strong  acceleration  prevents  the  transition  from 
proceeding.  The  beginning  of  the  rise  in  y  corresponds  to  K  drop¬ 
ping  below  3X1(F6.  In  low  FSTI  boundary  layers,  K> 3 
X  1CF6  leads  to  relaminarization  (Jones  and  Launder  [25]).  Fol¬ 
lowing  the  technique  of  Narasimha  [26],  as  modified  by  Volino 
and  Simon  [27],  the  function 

f(yPk)=(-  M 1  ~  yPk) ) 1/2  (8) 

can  be  computed  based  on  the  peak  intermittency  at  each  station 
and  plotted  versus  streamwise  location.  The  data  in  these  coordi¬ 
nates  tend  to  lie  along  a  straight  line.  The  line  may  be  extrapolated 
t°  /(7p*  =  0)  =  0  and  /(yp<r  =  0.99)  =  2.146,  corresponding  to  the 
beginning  and  end  of  transition  at  xs  =  0.29  m  and  xe  =  0.98  m, 
respectively.  Figure  5  shows  ypk  plotted  versus  dimensionless 
streamwise  location  within  transition.  Also  shown  is  the  theoreti¬ 
cal  curve 

yM  =  exp(-4.6(^)  )  (9) 

based  on  the  Dhawan  and  Narasimha  [28]  transition  model. 
Agreement  between  the  data  and  theory  is  good.  In  favorable 
pressure  gradients  some  differences  are  expected  and  observed  at 
low  intermittency,  in  a  region  Narasimha  [26]  referred  to  as  “sub- 
transition.” 

The  dimensionless  turbulent  spot  propagation  rate  may  be  com¬ 
puted,  following  the  development  of  Mayle  [1],  as 


4.6  v2U„ 
(xe~xs)2U3s 


(10) 


where  (7„  is  the  average  freestream  velocity  in  the  transition  re¬ 
gion.  In  the  present  case  ncr=  4.2X10~n.  In  agreement  with 
trends  reported  by  Mayle  [1]  for  other  favorable  pressure  gradient 
cases,  this  value  is  an  order  of  magnitude  lower  than  what  would 
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Fig.  8  Momentum  thickness  versus  streamwise  distance 


be  expected  for  a  zero-pressure  gradient  case  with  the  same  FSTI, 
indicating  a  longer  transition  zone  with  the  favorable  pressure 
gradient  than  with  a  zero-pressure  gradient. 

Mean  Velocity  Profiles.  Mean  velocity  profiles  for  the  ten 
measurement  stations  are  presented  in  Fig.  6  in  wall  coordinates. 
Figure  6(a)  shows  the  composite  (unconditioned)  profiles.  The 
profiles  rise  through  Station  6,  with  a  somewhat  laminar-like 
shape.  Downstream  of  Station  6,  the  profiles  assume  a  more 
turbulent-like  shape  and  by  Station  10,  where  transition  is  near 
completion  and  the  acceleration  rate  has  dropped  to  A"=l 
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Fig.  9  Shape  factor  versus  streamwise  distance 


Fig.  7  Mean  velocity  profile  for  Station  7,  ypk=  56% 


Fig.  10  Skin-friction  coefficient 
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X1CF6,  there  is  good  agreement  with  the  zero  pressure  gradient 
law  of  the  wall.  The  nonturbulent  zone  profiles  are  shown  in  Fig. 
6(b).  All  the  profiles  exhibit  a  laminar-like  shape,  even  at  the  end 
of  transition.  Figure  6(c)  shows  the  turbulent  zone  profiles.  All 
have  a  turbulent-like  shape,  and  the  last  five  stations  show  good 
agreement  with  the  zero-pressure  gradient  law  of  the  wall.  The 
high  FSTI  suppresses  the  wake  at  all  stations.  Figure  7  illustrates 
the  differences  between  the  composite,  nonturbulent,  and  turbu¬ 
lent  profiles  at  Station  7,  in  the  center  of  the  transition  zone. 
Distance  from  the  wall  is  normalized  on  the  composite  flow  899.5  ■ 
Velocities  in  the  turbulent  zone  are  clearly  higher  in  the  near  wall 
region  due  to  higher  levels  of  turbulent  mixing.  Also  shown  in 
Fig.  7  are  low  FSTI  calculations  for  laminar  and  fully  turbulent 


Fig.  11  Fluctuating  streamwise  velocity  profiles  in  wall  coor¬ 
dinates:  (a)  composite;  (ft)  nonturbulent;  (c)  turbulent 


Fig.  12  Fluctuating  streamwise  velocity  profile  for  Station  7, 

7Pk= 56% 


boundary  layers  with  the  same  pressure  gradient.  The  fully  turbu¬ 
lent  calculation  was  done  with  a  mixing  length  model.  The  turbu¬ 
lent  and  nonturbulent  zone  profiles  at  all  stations  differ  from  the 
low  FSTI  calculations.  Higher  levels  of  transport  lead  to  higher 
near  wall  mean  velocities  in  the  data  than  in  the  calculations.  This 
is  particularly  true  in  the  nonturbulent  zone.  The  increased  trans¬ 
port  in  the  nonturbulent  zone  makes  the  differences  between  the 
turbulent  and  nonturbulent  profiles  less  pronounced  in  the  data 
than  in  the  low  FSTI  calculations.  At  no  point  does  the  boundary 
layer  behave  as  if  it  were  laminar.  Differences  from  laminar  be¬ 
havior  have  been  reported  in  low  FSTI  transitional  flows  (e.g.,  [5], 
[15])  and  attributed  to  the  effect  of  the  turbulent  zone  on  the 
nonturbulent.  In  the  present  study,  however,  deviation  from  lami¬ 
nar  behavior  is  seen  even  at  Station  1,  where  the  boundary  layer  is 
nearly  all  nonturbulent.  The  deviation  from  laminar  behavior  must 
be  due  to  the  high  FSTI. 

Boundary  Layer  Growth.  Figure  8  shows  the  momentum 
thickness  as  a  function  of  streamwise  position.  The  composite 
boundary  layer  does  not  grow  through  the  first  five  stations  due  to 
the  strong  acceleration.  Momentum  thickness  increases  at  the 
downstream  stations  as  the  acceleration  weakens  and  transition 
proceeds.  In  the  nonturbulent  zone,  momentum  thickness  remains 
constant  at  all  stations.  The  turbulent  zone  momentum  thickness 
increases  continuously,  possibly  due  to  turbulent  entrainment  at 
the  edge  of  the  boundary  layer. 

The  shape  factor,  shown  in  Fig.  9,  is  an  indicator  of  the  state  of 
the  boundary  layer  with  respect  to  transition.  In  the  nonturbulent 
zone  it  drops  only  slightly  from  2.0  to  about  1.8.  A  low  FSTI 
laminar  boundary  layer  with  the  same  pressure  gradient  would 
have  a  shape  factor  of  2.4.  As  shown  in  Fig.  7,  the  high  FSTI 
enhances  mixing  and  makes  the  nonturbulent  zone  appear  less 
laminar-like.  In  the  turbulent  zone,  H  drops  from  about  1.5  to 
1.35.  The  low  FSTI  fully  turbulent  calculation  mentioned  above 
resulted  in  H  values  about  10%  higher  than  the  experimental  data. 
This  suggests  that  the  high  FSTI  promotes  greater  momentum 
transport  in  the  turbulent  zone,  but  the  effect  is  not  as  great  as  in 
the  nonturbulent  zone. 

Skin  Friction  Coefficients.  Skin  friction  coefficients  were 
computed  from  the  mean  velocity  profiles  and  are  shown  in  Fig. 
10.  Also  shown  for  reference  are  the  results  from  the  low  FSTI 
calculations.  The  skin  friction  coefficient  is  as  much  as  70% 
higher  in  the  turbulent  zone  than  in  the  nonturbulent  zone.  The 
composite  C'j  does  not  change  much  during  transition,  which  if 
viewed  alone,  might  suggest  that  there  is  little  difference  between 
the  nonturbulent  and  turbulent  zones.  When  viewed  with  the  con¬ 
ditional  sampling  results,  however,  it  is  clear  that  the  two  zones 
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are  quite  different.  In  both  zones  the  freestream  turbulence  results 
in  Cf  as  much  as  40%  higher  than  the  corresponding  low  FSTI 
calculation. 

Fluctuating  Velocity.  Figure  1 1  shows  u '  profiles  in  wall  co¬ 
ordinates.  The  composite  flow  data,  shown  in  Fig.  11(a),  are  typi¬ 
cal  of  transitional  and  turbulent  boundary  layers.  The  peak  in 
u' /uT  occurs  at  y+—  16,  and  the  magnitude  of  the  peak  is  between 
2  and  2.6.  The  peak  is  highest  at  Station  7,  in  the  middle  of 
transition,  due  in  part  to  the  unsteadiness  associated  with  the 
switching  between  turbulent  and  nonturbulent  states.  The  nontur- 
bulent  and  turbulent  zone  data  are  shown  in  Figs.  11(&)  and  11(c), 
respectively.  In  the  turbulent  zone  the  data  from  all  stations  col¬ 


lapse,  showing  self-similarity  throughout  the  transition  region. 
Figure  12  shows  the  u'  profiles  at  Station  7,  midway  through 
transition.  The  peak  in  u'  is  closer  to  the  wall  in  the  turbulent 
zone,  but  the  magnitudes  of  the  peaks  differ  by  only  about  13% 
between  the  two  zones. 

Figure  13  shows  v'  profiles.  The  composite  flow  data  are 
shown  in  Fig.  13(a).  As  expected  for  a  high  FSTI  boundary  layer, 
v '  drops  from  a  peak  in  the  near  wall  region  to  a  minimum,  and 
then  rises  to  the  freestream  value.  The  minimum  is  indicative  of  a 
damping  of  the  freestream  effect  by  the  wall.  Near  the  wall,  there 
is  not  much  change  in  v ' I ur  for  the  first  three  stations.  As  transi¬ 
tion  begins,  v' tur  rises  at  Stations  4  and  5.  Between  Stations  6 
and  10  there  is  little  change.  In  comparison  to  the  u'  profiles  of 
Fig.  11,  v '  is  lower  in  magnitude  and  shows  more  change  through 
the  transition  region.  Far  from  the  wall,  normalized  v'  drops  in 
the  streamwise  direction  due  to  the  increasing  value  of  uT.  The 
dimensional  value  of  v '  in  the  freestream  remains  essentially  con¬ 
stant  at  0.4  m/s  at  all  stations.  The  nonturbulent  and  turbulent  zone 
data  are  shown  in  Figs.  13 (b)  and  13(c).  As  with  u' ,  the  turbulent 
zone  profiles  collapse.  Figure  14  shows  the  v'  profiles  at  Station 
7.  In  contrast  to  the  u'  profiles  of  Fig.  12,  which  showed  similar 
magnitude  in  the  turbulent  and  nonturbulent  zones,  the  v '  magni¬ 
tude  is  nearly  twice  as  high  in  the  turbulent  zone  as  in  the  non¬ 
turbulent.  Volino  [17]  found  that  much  of  the  unsteadiness  in  u  is 
low-frequency  unsteadiness  induced  by  the  freestream  and  com¬ 
mon  to  both  zones.  When  freestream  eddies  buffet  the  boundary 
layer  and  push  higher  speed  fluid  toward  the  wall,  the  effect  is  an 
increase  in  u,  particularly  near  the  wall  where  dU/dy  is  large. 
This  effect  is  common  to  both  the  nonturbulent  and  turbulent 
zones  and  is  not  dependent  on  turbulence  produced  near  the  wall. 
It  explains  the  similarity  between  the  turbulent  and  nonturbulent 
u '  behavior.  There  is  no  similar  effect  on  v ' ,  so  the  u '  fluctua¬ 
tions  are  more  closely  related  to  turbulence  and  eddy  transport  in 
the  boundary  layer,  and  greater  differences  exist  between  the  tur¬ 
bulent  and  nonturbulent  v ' . 

Turbulent  Shear  Stress.  Profiles  of  the  turbulent  shear  stress 
are  shown  in  Fig.  15  in  wall  coordinates.  Similar  to  v ' luT  in  Fig. 
13(a),  there  is  little  change  in  the  composite  profiles  in  Fig.  15(a) 
for  the  first  three  stations.  As  transition  begins,  the  dimensionless 
shear  stress  rises.  In  the  nonturbulent  zone  (Fig.  15(b))  —u'v' 
values  are  low.  In  the  turbulent  zone,  Fig.  15(c)  shows  that  the 
data  from  beginning  to  end  of  transition  all  collapse  onto  a  single 
line.  To  compare  the  turbulent  and  nonturbulent  —u'v'  directly. 
Fig.  16  shows  the  Station  7  profiles.  In  contrast  to  u ' ,  which  had 
similar  magnitude  in  the  two  zones,  —u'v'  is  much  higher  in  the 


Fig.  13  Fluctuating  wall-normal  velocity  profiles  in  wall  coor-  Fig.  14  Fluctuating  wall-normal  velocity  profile  for  Station  7, 
dinates:  (a)  composite;  ( b )  nonturbulent;  (c)  turbulent  ypk=56% 
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Fig.  15  Reynolds  shear  stress  profiles  in  wall  coordinates:  (a)  Fig  17  Eddy  viscosity  profiles:  (a)  composite;  (to)  nonturbu- 
composite;  (to)  nonturbulent;  (c)  turbulent  lent.  turbulent 


turbulent  zone.  The  nonturbulent  —u'v'  is  not  zero,  however, 
indicating  some  eddy  transport  of  momentum  even  when  the 
boundary  layer  is  non-turbulent.  This  may  explain  the  deviation  of 
the  nonturbulent  mean  velocity  profiles  from  laminar  predictions 
and  the  enhancement  of  the  skin  friction  above  laminar  values 
shown  in  Figs.  7  and  10. 

The  collapse  of  the  turbulent  zone  data  in  Fig.  15(c),  also  seen 
in  u'  and  v'  in  Figs.  11(c)  and  13(c),  indicates  considerable  simi- 
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Fig.  16  Reynolds  shear  stress  profile  for  Station  7,  ypk= 56% 


larity  within  the  turbulent  zone.  The  turbulence  contained  within 
the  turbulent  spots  in  the  upstream  part  of  transition  (where  the 
intermittency  is  low  and  the  spots  occupy  only  a  small  fraction  of 
the  total  flow)  appears  to  be  very  similar  to  the  turbulence  in  the 
nearly  fully  turbulent  region  downstream.  This  similarity  may 
simplify  modeling  of  the  turbulent  zone.  Although  the  collapse  of 
the  data  in  the  coordinates  of  Figs,  ll(to),  13 (to),  and  15(to)  is  not 
quite  so  good,  similar  arguments  can  be  made  concerning  the 
self-similarity  of  the  nonturbulent  zone. 

Profiles  of  the  eddy  viscosity  are  shown  in  Fig.  17.  The  com¬ 
posite  profiles  in  Fig.  17(a)  show  that  em  increases  in  the  stream- 
wise  direction  as  the  transition  proceeds.  Comparison  of  the  non¬ 
turbulent  and  turbulent  zone  profiles  in  Figs.  11(b)  and  17(c) 
shows  that  the  eddy  viscosity  is  three  to  four  times  higher  in  the 
turbulent  zone.  While  the  difference  between  the  two  zones  is 
clear,  Fig.  11(b)  again  shows  that  there  is  significant  eddy  trans¬ 
port  in  the  nonturbulent  zone. 


Conclusions 

Conditional  sampling  was  successfully  performed  on  experi¬ 
mental  data  from  a  transitional  boundary  layer  subject  to  high 
freestream  turbulence  and  strong  acceleration.  Intermittency  func- 
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tions  based  on  the  instantaneous  streamwise  velocity  and  the  in¬ 
stantaneous  turbulent  shear  stress  agreed  well  and  produced  es¬ 
sentially  equal  conditional  sampling  results. 

Mean  velocity  profiles  differed  significantly  between  the  turbu¬ 
lent  and  nonturbulent  zones,  and  skin  friction  coefficients  in  the 
turbulent  zone  were  as  much  as  70%  higher  than  in  the  nonturbu¬ 
lent  zone.  The  u '  fluctuation  levels  did  not  differ  greatly  between 
the  turbulent  and  nonturbulent  zones,  but  v'  and  the  turbulent 
shear  stress  were  significantly  higher  in  the  turbulent  zone.  Within 
each  zone,  considerable  self-similarity  was  observed  in  all  turbu¬ 
lence  quantities  from  beginning  to  end  of  transition,  particularly  in 
the  turbulent  zone.  The  differences  between  the  two  zones  and  the 
similarity  within  each  zone  suggest  the  importance  of  properly 
modeling  the  transition  process  in  boundary  layer  prediction  and 
support  arguments  for  the  development  of  intermittency  based 
models. 

Although  the  turbulent  shear  stress  was  higher  in  the  turbulent 
zone,  nonturbulent  zone  eddy  transport  was  still  significant. 
In  both  zones  there  was  significant  deviation  from  low  FSTI 
predictions. 
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a  =  turbulent  spot  propagation  parameter 
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transition  end 
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